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Abstract 

The exclusive electroproduction of two pions in the mass range 
0.4 < M^^ < 2.5 GeV has been studied with the ZEUS detector at HERA 
using an integrated luminosity of 82 pb~^. The analysis was carried out in the 
kinematic range of 2 < < go GeV^ ?,2<W < 180 GeV and \t\ < 0.6 GeV^ 
where is the photon virtuality, W is the photon-proton centre-of-mass energy 
and t is the squared four-momentum transfer at the proton vertex. The two-pion 
invariant-mass distribution is interpreted in terms of the pion electromagnetic 
form factor, \F{Mt^tj)\, assuming that the studied mass range includes the 
contributions of the p, p' and p" vector-meson states. The masses and widths of 
the resonances were obtained and the dependence of the cross-section ratios 
a{p' — 7- mi)/a{p) and a{p" — t- n'n-)/a{p) was extracted. The pion form factor 
obtained in the present analysis is compared to that obtained in e~^e~ — t- tt+tt^. 



The ZEUS Collaboration 



H. Abramowicz^^'"'^, I. Abt'^^, L. Adamczyk^'^, M. Adamus^'^, R. AggarwaF''^, S. Antonelli^, 
P. Antonioli^, A. Antonov^^, M. Ameodo^°, D. Ashery^^, V. Aushev^'^'^'^'^, Y. Aushev,^'^'^'"", 
O. Bachynska^^, A. Bamberger^^, A.N. Barakbaev^^, G. Barbagli^^, G. Bari'^, F. Barreiro^^, 
N. Bartosik27'"^ D. Bartsch^, M. Basile"^, O. Behnkel^ J. Behv^^, U. Behrensl^ L. Bellagamba^, 
A. Bertolin^s, S. Bhadra^^, M. Bindi"^, C. Blohml^ V. Bokhonov^^.^, T. Bold^^^ k. Bondarenko^^, 

E. G. Boos^^, K. Borras^^, D. Boscherini^, D. Bot^^, I. Brock^, E. Brownson^^, R. Brugnera^", 
N. Briimmer'^'^, A. Bruni'^, G. Bruni^, B. Brzozowska^^ , P.J. Bussey^°, B. Bylsma'^^, A. Caldwell"^^, 
M. Capua^ R. Carlin^°, CD. Catterall5^ S. Chekanov\ J. Chwastowski^^'^, J. Ciborowski^^.ai^ 
R. Ciesielskii^'f , L. Cifarelli^ F. Cindolo^, A. Contin^ A.M. Cooper-Sarkar^s, N. Coppola^^''', 
M. Corradi^, F. Corriveau^i, M. Costa^^ G. D'Agostini^^ F. Dal Corso^^, J. del Peso^o, 
R.K. Dementiev^^, S. De Pasquale^'"' , M. Demck\ R.C.E. Devenish^'^, D. Doburl9'^ B.A. Dolgo- 
shein G. Dolinska^'^'^'^, A.T. Doyle^o, V. Drugakov^^ L.S. Durkin^^, S. Dusini^s, Y. Eisenberg5^ 
P.F. Ermolov ^^'^ A. Eskreys ^^'t, S. Fang^^'*, S. Fazio^, J. Ferrando^^, M.I. Ferrero^^, J. Figiel^^, 
M. Forrest^o.'', B. Foster^s.'^'^, G. Gachi^, A. Galas^^, E. Gallo^^ A. Garfagnini^o , A. Geiserl^ 

I. Gialas^i''", L.K. Gladilin^^.^^, D. Gladkov^^, C. Glasman^o, 0. Gogota26.27^ Yu.A. Golubkov^^, 
P. Gottlicher^^'-' , I. Grabowska-Bold^'^, J. Grebenyuk^^, I. Gregor^^, G. Grigorescu^^ , G. Grzelak^^, 
O. Gueta^^ E. Gurvich^^ M. Guzik^^, C. GwenW^''^^ T. Haasl^ W. Hainl^ R. Hamatsu^^ 
J.C. Hart^"^, H. Hartmann^ G. Hartner^^, E. Hilger^ D. Hochman^s, R. HorF, K. Horton^S''^/, 
A. Hiittmanni^ Z.A. Ibrahimi°, Y. Iga^^ R. Ingbir'^^ M. Ishitsuka^^, H.-P. Jakob^ F. Januschek^^ 
T.W. Jones^^ M. Jungst^ I. Kadenko^^, B. Kahle^^ S. Kananov^^ T. Kanno'^^ U. Karsllon5^ 

F. Karstens^^'*, I.I. Katkov^^-*^, M. Kaur^ P. Kaur^'^ A. Keramidas^e, L.A. Khein^^, J.Y. Kim^ 
D. Kisielewska^'^ , S. Kitamura^®'"-' , R. Klanner^^, U. Klein^^'', E. Koffeman^^, P. Kooijman^^, 
le. i.A. Korzhavina^^'^^, A. Kotanski^^'-^, U. Kotz^s, H. Kowalski^s, Q. Kuprash^s, 
M. Kuze^^ A. Lee^^ B.B. Levchenko3^ A. Levy^^ V. Libovl^ S. Limentani'^o , T.Y. Ling^^, 
M. Lisovyi^^, E. Lobodzinska^^, W. Lohmann^^, B. Lohr^^, E. Lohrmann^^, K.R. Long^^, 
A. Longhin^s, D. Lontkovskyi^^, O.Yu. Lukina^^, J. Maeda^'^'''*, S. Magill\ I. Makarenko^^, 
J. Malka^^, R. Mankel^^, A. Margotti^, G. Marini^^^ Martin^^ A. Mastroberardino^, 
M.C.K. Mattingly2, I.-A. Melzer-Pellmann^^, S. Mergelmeyer^ , S. Miglioranzi^^'"^, F. Mohamad 
Idris^°, V. Monaco"^^, A. Montanari^^, J.D. Moms^'^ K. Mujkic^^-'", B. Musgrave\ K. Nagano^"^, 
T. Namsoo^^'^ R. Nania^, A. Nigro^^^ Y. Ning", T. Nobe^^ U. Noor^^ D. Notzl^ R.J. Nowak^^^ 
A.E. Nuncio-Quiroz^, B.Y. Oh'^i, N. Okazaki^^ K. Oliver^s, K. Olkiewiczi^, Yu. Onishchuk^^, 

K. Papageorgiu^i, A. Parenti^^, E. Paul^ J.M. Pawlak^^^ g Pawlik^^ p. G. Pelfer^s, A. Pellegrino^*' , 
W. Perlanski^^'""", H. Perrey^^ K. Piotrzkowski^^ , P. Plucihski^^'''", N.S. Pokrovskiy^^ , A. Polini^, 
A.S. Proskuryakov^^, M. Przybycien^^, A. Raval^^, D.D. Reeder^^, B. Reisert^^, Z. Ren^\ 
J. Repond^, Y.D. Ri^^.afc^ Robertson^s, P. Roloffi^'™, I. Rubinskyi^, M. Ruspa^", R. Sacchi^^^ 
A. Salii27, U. Samson^, G. Sartorelli^, A.A. Savm^e, D.H. Saxon^o, M. Schioppa^, S. Schlenstedt^^, 
P. Schleper^^, W.B. Schmidke^^, U. Schneekloth^^ , V. Schonberg^, T. Schorner-Sadenius^^, 
J. Schwartz^i, F. Sciulli", L.M. Shcheglova^^ , R. Shellzadi^ S. Shimizu^^'"*, I. Singh^•^ I.O. Skillicorn' 
W. Slominskil^ W.H. Smith5^ V. Sola'^^ A. Solano'^^ D. Son^^, V. Sosnovtsev^s, A. Spiridonoyi^'P, 
H. Stadie22, L. Stanco^s, A. Stern^^^ T.P. Stewart^i, A. Stifutkin^^, p. Stopa^^, S. Suchkov^s, 



I 



G. Susinno®, L. Suszycki^^, J. Sztuk-Dambietz^^ , D. Szuba^^, J. Szuba^^-'', A.D. Tapper^^, 
E. Tassi^''^, J. Terron^o, T. Theedt^^ H. Tiecke3^ K. Tokushuku^^''' , O. Tomalak^'^, J. Tomaszewska^^''' , 
T. Tsurugai32, M. Turcato^^, T. Tymieniecka^'^'"'" , M. Vazquez^^.™, A. Verbytskyil^ O. Viazlo^^.s?^ 
N.N. Vlasov^S'", O. Volynets^^, R. Walczak^s, W.A.T. Wan Abdullahio, J.J. Whitmore'^i''^^', 
L. Wiggers^e, M. Wing^^^ M. Wlasenko^ G. Wolf^^ H. Wolfe^^ K. Wrona^^ A.G. Yagiies- 
Molina-'^^, S. Yamada^^, Y. Yamazaki^"^'^ , R. Yoshida-^, C. Youngman-'^^ , A.F. Zarnecki^'^, L. Zawiejski-*^^, 
O. Zenaiev^s, W. Zeuner^S'"^, B.O. Zhautykov^^, N. Zhmak^e.^, C. Zhou3^ A. Zichichi'^, Z. Zolkapli^", 
M. Zolko27, D.S. Zotkin34 



II 



Argonne National Laboratory, Argonne, Illinois 60439-4815, USA ^ 

Andrews University, Berrien Springs, Michigan 49104-0380, USA 

INFN Bologna, Bologna, Italy ^ 

University and INFN Bologna, Bologna, Italy ^ 

Physikalisches Institut der Universitdt Bonn, Bonn, Germany ^ 

H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom ^ 

Panjab University, Department of Physics, Chandigarh, India 

Calabria University, Physics Department and INFN, Cosenza, Italy ^ 

Institute for Universe and Elementary Particles, Chonnam National University, 

Kwangju, South Korea 

Jabatan Fizik, Universiti Malaya, 50603 Kuala Lumpur, Malaysia ^ 

Nevis Laboratories, Columbia University, Irvington on Hudson, New York 10027, 

USA ^ 

The Henryk Niewodniczanski Institute of Nuclear Physics, Polish Academy of 
Sciences, Krakow, Poland 

ACH-University of Science and Technology, Faculty of Physics and Applied Com- 
puter Science, Krakow, Poland ^ 

Department of Physics, Jagellonian University, Cracow, Poland 
Deutsches Elektronen- Synchrotron DESY, Hamburg, Germany 
Deutsches Elektronen- Synchrotron DESY, Zeuthen, Germany 
INFN Florence, Florence, Italy ^ 
University and INFN Florence, Florence, Italy ^ 

Fakultdt fiir Physik der Universitdt Freiburg i.Br., Freiburg i.Br., Germany 
School of Physics and Astronomy, University of Glasgow, Glasgow, United King- 
dom ^ 

Department of Engineering in Management and Finance, Univ. of the Aegean, 
Chios, Greece 

Hamburg University, Institute of Experimental Physics, Hamburg, Germany ^ 
Imperial College London, High Energy Nuclear Physics Group, London, United 
Kingdom ^ 

Institute of Particle and Nuclear Studies, KEK, Tsukuba, Japan 
Institute of Physics and Technology of Ministry of Education and Science of Kaza- 
khstan, Almaty, Kazakhstan 

Institute for Nuclear Research, National Academy of Sciences, Kyiv, Ukraine 
Department of Nuclear Physics, National Taras Shevchenko University of Kyiv, 
Kyiv, Ukraine 

Kyungpook National University, Center for High Energy Physics, Daegu, South Ko- 
rea ^ 

Institut de Physique Nucleaire, Universite Catholique de Louvain, Louvain-la-Neuve, 
Belgium ^ 

Departamento de Fisica Teorica, Universidad Autdnoma de Madrid, Madrid, 
Spain 

Department of Physics, McGill University, Montreal, Quebec, Canada H3A 2T8 ^ 
Meiji Gakuin University, Faculty of General Education, Yokohama, Japan ^ 



III 



Moscow Engineering Physics Institute, Moscow, Russia ^ 

Moscow State University, Institute of Nuclear Physics, Moscow, Russia ^ 

Max- Planck- Institut fiir Physik, Miinchen, Germany 

NIKHEF and University of Amsterdam, Amsterdam, Netherlands 

Physics Department, Ohio State University, Columbus, Ohio 43210, USA ^ 

Department of Physics, University of Oxford, Oxford, United Kingdom ^ 

INFN Padova, Padova, Italy ^ 

Dipartimento di Fisica dell' Universitd and INFN, Padova, Italy ^ 

Department of Physics, Pennsylvania State University, University Park, 

Pennsylvania 16802, USA ^ 

Polytechnic University, Sagamihara, Japan 

Dipartimento di Fisica, Universitd 'La Sapienza' and INFN, Rome, Italy ^ 
Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, United Kingdom ^ 
Raymond and Beverly Sackler Faculty of Exact Sciences, School of Physics, 
Tel Aviv University, Tel Aviv, Israel ^ 

Department of Physics, Tokyo Institute of Technology, Tokyo, Japan ^ 
Department of Physics, University of Tokyo, Tokyo, Japan 
Tokyo Metropolitan University, Department of Physics, Tokyo, Japan ^ 
Universitd di Torino and INFN, Torino, Italy ^ 

Universitd del Piemonte Orientale, Novara, and INFN, Torino, Italy ^ 
Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 
1A7^ 

Physics and Astronomy Department, University College London, London, United 
Kingdom ^ 

Faculty of Physics, University of Warsaw, Warsaw, Poland 
National Centre for Nuclear Research, Warsaw, Poland 

Department of Particle Physics and Astrophysics, Weizmann Institute, Rehovot, 
Israel 

Department of Physics, University of Wisconsin, Madison, Wisconsin 53706, USA ^ 
Department of Physics, York University, Ontario, Canada M3J 1P3 ^ 



IV 



supported by the US Department of Energy 

^ supported by the ItaHan National Institute for Nuclear Physics (INFN) 

supported by the German Federal Ministry for Education and Research (BMBF), 

under contract No. 05 H09PDF 
^ supported by the Science and Technology Facilities Council, UK 

^ supported by an FRGS grant from the Malaysian government 

^ supported by the US National Science Foundation. Any opinion, findings and con- 
clusions or recommendations expressed in this material are those of the authors and 

do not necessarily reflect the views of the National Science Foundation. 
^ supported by the Polish Ministry of Science and Higher Education as a scientific 

project No. DPN/N188/DESY/2009 
^ supported by the Polish Ministry of Science and Higher Education and its grants 

for Scientific Research 
^ supported by the German Federal Ministry for Education and Research (BMBF), 

under contract No. 05h09GUF, and the SFB 676 of the Deutsche Forschungsge- 

meinschaft (DFG) 

supported by the Japanese Ministry of Education, Culture, Sports, Science and 
Technology (MEXT) and its grants for Scientific Research 
^ supported by the Korean Ministry of Education and Korea Science and Engineering 
Foundation 

^ supported by FNRS and its associated funds (HSN and FRIA) and by an Inter- 
University Attraction Poles Programme subsidised by the Belgian Federal Science 
Policy Office 

^ supported by the Spanish Ministry of Education and Science through funds provided 
by CICYT 

^ supported by the Natural Sciences and Engineering Research Council of Canada 
(NSERC) 

'-^ partially supported by the German Federal Ministry for Education and Research 
(BMBF) 

^ supported by RF Presidential grant N 4142.2010.2 for Leading Scientific Schools, 
by the Russian Ministry of Education and Science through its grant for Scientific 
Research on High Energy Physics and under contract No. 02. 740. 11. 0244 

^ supported by the Netherlands Foundation for Research on Matter (FOM) 

^ supported by the Israel Science Foundation 



V 



" now at University of Salerno, Italy 

^ now at Queen Mary University of London, United Kingdom 

^ also funded by Max Planck Institute for Physics, Munich, Germany 

also Senior Alexander von Humboldt Research Fellow at Hamburg University, Insti- 
tute of Experimental Physics, Hamburg, Germany 
^ also at Cracow University of Technology, Faculty of Physics, Mathemathics and 

Applied Computer Science, Poland 
f supported by the research grant No. 1 P03B 04529 (2005-2008) 
^ now at Rockefeller University, New York, NY 10065, USA 
^ now at DESY group FS-CFEL-1 

* now at Institute of High Energy Physics, Beijing, China 
^ now at DESY group FEB, Hamburg, Germany 

^ also at Moscow State University, Russia 

' now at University of Liverpool, United Kingdom 

now at CERN, Geneva, Switzerland 
" also affiliated with Universtiy College London, UK 
° now at Goldman Sachs, London, UK 

^ also at Institute of Theoretical and Experimental Physics, Moscow, Russia 

<? also at FPACS, AGH-UST, Cracow, Poland 

^ partially supported by Warsaw University, Poland 

* now at Istituto Nucleare di Fisica Nazionale (INFN) , Pisa, Italy 

* now at Haase Energie Technik AG, Neumiinster, Germany 
now at Department of Physics, University of Bonn, Germany 

" now at Biodiversitat und Klimaforschungszentrum (BiK-F), Frankfurt, Germany 

also affiliated with DESY, Germany 
^ also at University of Tokyo, Japan 
y now at Kobe University, Japan 
^ supported by DESY, Germany 
^ deceased 

aa jj^ember of National Technical University of Ukraine, Kyiv Polytechnic Institute, 
Kyiv, Ukraine 

""^ member of National University of Kyiv - Mohyla Academy, Kyiv, Ukraine 
"•'^ partly supported by the Russian Foundation for Basic Research, grant 11-02-91345- 
DFG_a 

"'^ Alexander von Humboldt Professor; also at DESY and University of Oxford 

STFC Advanced Fellow 
"■^ nee Korcsak-Gorzo 

This material was based on work supported by the National Science Foundation, 

while working at the Foundation. 
"'^ also at Max Planck Institute for Physics, Munich, Germany, External Scientific 

Member 

now at Tokyo Metropolitan University, Japan 
"■^ now at Nihon Institute of Medical Science, Japan 



VI 



now at Osaka University, Osaka, Japan 
also at Lodz University, Poland 
member of Lodz University, Poland 

now at Department of Physics, Stockholm University, Stockholm, Sweden 
also at Cardinal Stefan Wyszyhski University, Warsaw, Poland 



VII 



1 Introduction 



Exclusive electroproduction of vector mesons takes place through a virtual photon 7* 
by means of the process 7*p — )■ Vp. At large values of the centre-of-mass energy, W, 
this is usually viewed as a three-step process; the virtual photon 7* fluctuates into a qq 
pair which interacts with the proton through a two-gluon ladder and hadronizes into a 
vector meson, V. The production of ground-state vector mesons, V = p,u,(f>, J/ip,T, 
which are IS" triplet qq states, has been extensively studied at HERA, particularly in 
several recent publications [IHZ])- As the virtuality, Q^, of the photon increases, the 
process becomes hard and can be calculated in perturbative Quantum Chromodynamics 
(pQCD). Furthermore, by varying Q^, and thus the size of the qq pair, sensitivity to 
the vector-meson wave-function can be obtained by scanning it at different qq distance 
scales. Expectations in the QCD framework vary from calculations based only on the mass 
properties and typical size of the qq inside the vector- meson [8l|9], to those based on the 
details of the vector meson wave-function dependence on the size of the qq pair [TUHTi] . 
The approaches differ in their predictions for the dependence of the cross sections for 
excited vector-meson states and their ratio to their ground state. 

The only radially excited 25* triplet qq state studied at HERA so far has been the ip{'2S) 
state [15]. In this study, only the photoproduction reaction was investigated and the 
low cross-section ratio of ipi'iS) to the ground-state J/ip supported the existence of a 
suppression effect, expected if a node in the ip{2S) wave-function is present. 

Other excited vector-meson states, in particular those consisting of light quarks, can be 
used to study the effect caused by changing the scanning size. Exclusive n~^n~ production 
has been measured previously in the annihilation process e^e~ — )■ vr+vr" [16], as well as in 
photoproduction [T7]. The tt+tt" mass distribution shows a complex structure in the mass 
range 1-2 GeV. Evidence for two excited vector- meson states has been established [T8|[T9] ; 
the p'(1450) is assumed to be predominantly a radially excited 25* state and the p"(1700) 
is an orbitally excited 2D state, with some mixture of the S and D waves [20]. In addition 
there is also the p3(1690) spin-3 meson [21] which has a tttt decay mode. The two-pion 
decay mode of these resonances is related [2^125] to the pion electromagnetic form factor. 

In this paper, a study of exclusive electroproduction of two pions, 

■y*P n^n^p, (1) 

is presented in the two-pion mass range 0.4 < M^rTr < 2.5 GeV, in the kinematic range 
2 < g2 < 80 GeV^, 32 < W < 180 GeV and \t\ < 0.6 GeV^ where t is the squared 
four-momentum transfer at the proton vertex. The Mt^t^ system consists of a resonance 
part and a non-resonant background. The resonances are described by the pion form 



1 



factor. The contributions of the three vector-mesons p, p' and p" are extracted and their 
relative rates as a function of are discussed in terms of QCD expectations. 



2 The pion form factor 



The two-pion invariant-mass distribution of Eq. [H after subtraction of the non-resonant 
backgrouncQ, can be related to the pion electromagnetic form factor, Ft^{Mt^t^\ through 
the following relation 



(2) 



There are several parameterizations of the pion form factor usually used for fitting the 
TT+vr" mass distribution; the Kuhn-Santamaria (KS) [21], the Gounaris-Sakurai (GS) [25] 
and the Hidden Local Symmetry (HLS) [26l|27| parameterizations. In this paper, results 
based on the KS parameterization are presented. 

In the mass range M^r^ < 2.5 GeV, the KS parameterization of the pion form factor 
includes contributions from the p(770), p'(1450) and p"(1700) resonance^. 



(3) 



1 + /3 + 7 

Here /3 and 7 are relative amplitudes and BWy is the Breit-Wigner distribution which 
has the form 



BWviM^ 



(4) 



M^-M^^-tMvTv{M^.y 

where My and Tv{Mt^t,) are the vector-meson mass and momentum- dependent width, 
respectively. The latter has the form 



pJMr 



V) 



Ml 
ML 



(5) 



where Ty is the width of the V meson at M^^ = My, p^M^^ = 1/2 a/M^^ - 4M^ is the 
pion momentum in the tt+tt^ centre-of-mass frame, p^(Mv') is the pion momentum in the 
l^-meson rest frame, and M^^ is the pion mass. 



^ This is assumed not to interfere with the resonance signal. 

^ This analysis cannot distinguish between p3(1690) and p"(1700). Theoretical calculations estimate 
the contribution of p3(1690) to be either one order of magnitude [9] or 2-5 times [28] smaller than that 
of the p"(1700). 
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3 Experimental set-up 



The analyzed data were collected with the ZEUS detector at the HERA collider in the 
years 1998-2000, when 920 GeV protons collided with 27.5 GeV electrons or positrons. 
The sample used for this study corresponds to 81.7 pb^^ of which 65.0 pb~^ were collected 
with an e"*" and the rest with an e~ beanif|. 

A detailed description of the ZEUS detector can be found elsewhere [29] . A brief outline 
of the components that are most relevant for this analysis is given below. 

The charged particles were tracked in the central tracking detector (CTD) [30] which 
operated in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The 
CTD consisted of 72 cylindrical drift chamber layers, organized in nine superlayers cov- 
ering the polar-angl^ region 15° < 9 < 164°. The transverse-momentum resolution for 
full-length tracks was cr{pT)/pT = 0.0058pT © 0.0065 © 0.0014/pt, with pr in GeV. 

The scattered electron was identified in the high-resolution uranium-scintillator calorime- 
ter (GAL) [31] which covered 99.7% of the total solid angle and consisted of three parts: 
the forward (FGAL), the barrel (BGAL) and the rear (RGAL) calorimeters. Each part 
was subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion (EMG) and either one (in RGAL) or two (in BGAL and FGAL) hadronic sec- 
tions (HAG). The GAL energy resolution, as measured under test-beam conditions, was 
a{E)/E = 0.18/y/E for electrons and a{E)/E = 0.35//E for hadrons, with E in GeV. 

The position of the scattered electron was determined by combining information from the 
GAL, the small-angle rear tracking detector [32] and the hadron-electron separator [33] . 

The luminosity was measured from the rate of the bremsstrahlung process ep — 67 p. 
The photon was measured in a lead-scintillator calorimeter [MH36] placed in the HERA 
tunnel at Z = — 107 m. 

4 Data selection and reconstruction 

The online event selection required an electron candidate in the GAL, along with the 
detection of at least one and not more than six tracks in the GTD. 

In the offline selection, the following further requirements were imposed: 

^ From now on, electrons and positrons will be both referred to as electrons in this paper. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The polar angle, 0, is 
measured with respect to the proton beam direction. 
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• the presence of a scattered electron, with energy in the CAL greater that 10 GeV 
and with an impact point on the face of the RCAL outside a rectangular area of 
26.4 X 16 cm^ in the X — Y plane; 

• the Z coordinate of the interaction vertex was within ± 50 cm of the nominal inter- 
action point; 

• in addition to the scattered electron, the presence of exactly two oppositely charged 
tracks. Both tracks have to be associated with the reconstructed vertex, each having 
pseudorapidity |?7| less than 1.75 and transverse momentum greater that 150 MeV. 
This ensures high reconstruction efficiency and excellent momentum resolution in the 
CTD. These tracks were treated in the following analysis as a tt'^tt" pair; 

• — > 45 GeV, where E — Pz = Ylii-^i ~ ^Zi) and the summation is over the 
energy Ei and longitudinal momentum P^- of the final-state electron and pions. This 
cut excludes events with high-energy photons radiated in the initial state; 

• events with any energy deposit larger than 300 MeV in the CAL, not associated with 
the pion tracks (so-called 'unmatched islands'), were rejected. 

The following kinematic variables are used to describe the exclusive production of a tt+tt^ 
pair: 

• Q^, the four-momentum squared of the virtual photon; 

• W^, the squared centre-of-mass energy of the photon-proton system; 

• Mtttt; the invariant mass of the two pions; 

• t, the squared four- momentum transfer at the proton vertex; 

• the angle between the tt+tt^ production plane and the positron scattering plane 
in the 7*p centre-of-mass frame; 

• Oh and (j)hi the polar and azimuthal angles of the positively charged pion in the s- 
channel helicity frame [3Z] of the tt+tt". 

The kinematic variables were reconstructed using the so-called 'constrained' method |38j . 
which uses the momenta of the decay particles measured in the CTD and the reconstructed 
polar and azimuthal angles of the scattered electron. The analysis was restricted to 
the kinematic region 2 < < GeV^ ?>2 < W < 180 GeV, |t| < 0.6 GeV^ and 
0.4 < Mjr^ < 2.5 GeV. The lower mass range excludes reflections from the cf) — ?■ K^K~ 
decays and the upper limit excludes the J /ip /^^/^~) e"'"e~ decays with its radiative tail. 

The above selection yielded 63517 events for this analysis. 

The above cuts do not eliminate events in which the proton dissociates into a low-mass 
final state, the products of which disappear down the beam pipe. This contribution. 
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estimated [2] to be about 20% in the range of this analysis, was found to be and W 
independent. Its presence does not affect the conclusions of this analysis. 



5 Monte Carlo simulation 

The program Zeusvm [39] interfaced to Heracles4.4 [ID] was used. The effective dis- 
tributions of Q^, W and |t| were parameterized to reproduce the data. The mass and 
angular distributions were generated uniformly and the MC events were then iteratively 
reweighted using the results of the analysis. 

The generated events were passed through a full simulation of the ZEUS detector based 
on Geant 3.21 [H] and processed through the same chain of selection and reconstruction 
procedures as the data, accounting for trigger as well as detector acceptance and smearing 
effects. The number of simulated events after reconstruction was approximately seven 
times greater than the number of reconstructed data events. 

A detailed comparison between the data and the Zeusvm MC distributions for the mass 
range 0.65 < M^^r < 1-1 GeV has been presented elsewhere [2]. Some examples for the 
mass range 1.1 < M^^r < 2.1 GeV are shown here. The transverse momentum, p^-, of 
the 7r+ and the tt" particles for different ranges of and M^tt are presented in Fig. [1] 
Figure |2] shows the Q^, cos 9h,(ph, and $/i distributions for events selected within 

the mass ranges 1.1 < Mj^t^ < 1.6 GeV, while Fig. |3] shows those distributions for the 
mass range 1.6 < M^tt < 2.1 GeV. All measured distributions are well described by the 
MC simulations. 



6 The 7V7V mass fit 



The TT+vr" mass distribution, after acceptance correction determined from the above MC 
simulation, is shown in Fig. HI A clear peak is seen in the p mass range. A small shoulder 
is apparent around 1.3 GeV and a secondary peak at about 1.8 GeV. 

The two-pion invariant-mass distribution was fitted, using the least-square method 
as a sum of two terms, 



dN{M^ 



A 1 



AMI 
ML 



Mo 



(6) 



where A is an overall normalization constant. The second term is a parameterization of 
the non-resonant background, with constant parameters B, n and Mq = 1 GeV. The other 
parameters, the masses and widths of the three resonances and their relative contributions 
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/9 and 7, enter through the pion form factor, (Eq. |3]). The fit, which includes 11 
parameters, gives a good description of the data (x^/ndf=28. 8/24=1. 2). The result of 
the fit is shown in Fig. H] together with the contribution of each of the two terms of Eq. [61 
The p and the p" signals are clearly visible. The negative interference between all the 
resonances results in the p' signal appearing as a shoulder. To illustrate this better, the 
same data and fit are shown in Fig. [5] on a linear scale and limited to Mtttt > 1.2 GeV, 
with separate contributions from the background, the three resonant amplitudes as well 
as their total interference term. 

The fit parameters are listed in Table [U Also listed are the mass and width parameters 
from the Particle Data Group (PDG) |13]. The masses and widths of the p and the p" 
as well as the width of the p' agree with those listed in the PDG, while there is about 100 
MeV difference between the PDG value and the fitted mass of the p' . It should however 
be noted that the value quoted by PDG is an average over many measurements having a 
large spread (1265 ± 75 up to 1424 ± 25 MeV for the vrvr decay mode) in this mass range. 

The measured negative value of (3 and positive value of 7 implies that the relative signs 
of the amplitudes of the three resonances p, p' and p" are +,—,+, respectively. A similar 
pattern was observed in e^e~ — ?■ ti^ti" and r-decay experiments [15H53], which also 
showed a dip in the mass range around 1.6 GeV, resulting from destructive interference. 
There is a single experiment where a constructive interference was obtained around 1.6 
GeV, namely 7p — t- tx'^ti^p [T7], a result which is not understood [T2] . 

In the mass fits above it was assumed that the relative amplitudes /3 and 7 are real. In 
order to test this assumption, the fit was repeated allowing them to be complex. The 
pion form factor was re-written in the form 

1^, ^ BWp{M^^) + /3o ■ eM^^ 12) BW p. {M^^) + 70 ■ e^^{i^i^)BWpn{M^^) 

1 + PO + 7o 

where /3o and 70 are real numbers and two additional fit parameters, $12 and $13, are 
the corresponding phase shifts. The value of the phase-shifts obtained from the fit were 
$12 = 3.2 ± 0.2 rad and $13 = 0.1 ± 0.2 rad, supporting the assumption of the real nature 
of the relative amplitudes. 



7 Systematic uncertainties 

The systematic uncertainties of the fit parameters were evaluated by varying the selection 
cuts and the MC simulation parameters. Motivation for the variation in cuts used below 
can be found in a previous ZEUS analysis [2] . The following selection cuts were varied: 

• the E — Pz cut was changed within the resolution of ± 3 GeV; 
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• the pt threshold for the pion tracks (default 0.15 GeV) was increased to 0.2 GeV and 
the |?7| cut on the two pion tracks was changed (default 1.75) by ±0.25; 

• the required maximum distance of closest approach of the two extrapolated pion tracks 
to the matched island in the CAL was changed from 30 cm to 20 cm; 

• the Z- vertex cut was varied by ±10 cm; 

• the energy threshold for an unmatched island (elasticity cut) was changed by ±50 MeV; 

• the bin size in the fitted mass distribution (default 60 MeV) was varied by ± 20 MeV; 

• the mass range was narrowed to 0.5 < Af^r^ < 2.3 GeV; 

• the \t\ cut was varied by ±0.1 GeV^; 

• the W range was changed to 35 < < 190 GeV; 

• the cosOfi range was changed to | cos6'/i| < 0.9; 

• the dependence in the MC was varied by changing the Q^-dependent S value by 
± 0.03; 

• the exponential t distribution in the MC was reweighted by changing the nominal 
Q^-dependent slope parameter 6 by ± 0.5 GeV~^; 

• the exponent of the distribution parameterization in the MC was changed by ±0.05. 

The largest variations were observed for 7, T{p") and (3. The value of T{p") changes by 
7% when the elasticity cut is varied. The restriction of the phase space in the fitted mass 
range leads to a change of the value of /3 by —5.2% while for 7, restricting the | cos 9h\ 
range leads to a change of —8%. In addition, another form of background in Eq. [6|, with 
an added exponential term, was investigated. It gave a very similar result in the mass 
range of this analysis and therefore no additional uncertainty was assigned to the form of 
the fitted mass curve. 

All the systematic uncertainties were added in quadrature. The combined systematic 
uncertainties are included in Table [TJ 

8 Decay angular distributions 

Decay angular distributions can be used to determine the spin density-matrix elements 
of a resonance [371111]. In the present case we study three resonances, all in a = 1^ 
state. However, the decay angular distribution in a given mass bin is affected by the 
background contribution which does not necessarily have the same quantum numbers as 
the resonance. Given the above, only the distribution of the polar angle 6h, defined as 
the polar angle of the positively charged pion in the helicity frame, was studied. 
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The distribution of cos 9h is shown in Fig. [6] for different mass bins; its shape is clearly 
mass dependent. In order to study the mass dependence further, the angular distribution 
of the polar helicity angle, W{cos9h) was parameterized as 

W{coseh) (x[l-r + {?,r-l)cos'^eh], (8) 

and fitted to the data. The mass dependence of the resulting parameter r is shown 
in Fig. [71 In the mass range M^^^ < 1.1 GeV, r shows the dependence seen for the Tqq 
density matrix in the p region [2] . Indeed this region is dominated by exclusive production 
of p and therefore r = Tqq. In that case, r can be interpreted as ai/atot, assuming s- 
channel helicity conservation (SCHC). Here 0"^ is the cross section for producing p by a 
longitudinally polarized photon, and atot = ctl + ct, with ar the production cross section 
by transversely polarized photons. The results shown here for the p region are in excellent 
agreement with the values given in an earlier ZEUS paper [2]. 

The structure seen for M^^r > 1.1 GeV is not easy to interpret, however the dip observed 
around 1.3 GeV and the enhancement at 1.6 GeV seem to follow the location of the 
resonances determined from the mass distribution. 



9 dependence of the pion form factor 

The dependence of the relative amplitudes was determined by performing the fit to 
M^TT in three regions, 2-5, 5-10 and 10-80 GeV^. The masses and widths of the three 
resonances were fixed to the values found in the overall fit and listed in Table [U The 
results are shown in Fig. [HI A reasonable description of the data is achieved in all three 
regions. The corresponding values of /3 and 7 are given in Table [2l The absolute value 
of /3 increases with while the value of 7 is consistent with no dependence, within 
large uncertainties. 

Figure [HI shows the curves representing the pion form factor, |F^(M^^)p, as obtained in 
the present analysis for the three Q ranges: 2-5, 5-10, 10-80 GeV^ Also shown are 
results obtained in the time-like regime from the reaction e"'"e~ tt^tc^ . In general, the 
features of the |F^(M^7r)P distribution observed here are also observed in e+e", i.e., the 
prominent p peak, a shoulder around the p' and a dip followed by an enhancement in the 
p" region. Above the p region, where the interference between the p' and the p" starts to 
dominate, there is a dependence of |F7r(M,r7r)P on Q^, with the results from the lowest 
range closest to those from e~^e~ . However, in the region of the p peak, shown in Fig. [TOl 
the pion form-factor [^^^(M^Tr)^ is highest at the highest Q^, as in the p'-p" interference 
region, while the e^e~ data are higher than those in the highest range. They are equal 
within errors for ikL^ > 1.8 GeV. 
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10 Cross-section ratios as a function of 



The dependence of the p by itself is given elsewhere [2]. Since the tttt branching ratios 
of p' and p" are poorly known, the ratio Ry defined as 

a{V) ■ Br\V ^ tttt) 
= ' 

has been measured, where a is the cross section for vector-meson production and Br(y — )■ 
tttt) is the branching ratio of the vector meson V{p',p") into tttt. The ratio Ry niay be 
directly determined from the results of the M^^r niass fit, 

R, = f3'^-f V=7^^, (10) 

where 

/y= y dM^^\BWv{M^.)\^ (11) 

and the integration is carried out over the range 2M7r < M^tt < My + SFy. 

Figure ITT] shows and Table E] lists the ratio Ry for = p',p", as a function of Q^. Owing 
to the large uncertainties of Rp", no conclusion on its behaviour can be deduced, 
whereas Rp> clearly increases with Q^. This rise has been predicted by several models [H 
UnillSllSllES]- The suppression of the 25* state (p') is connected to a node effect which 
results in cancellations of contributions from different impact-parameter regions at lower 
Q^, while at higher the effect vanishes. 



11 Summary 

Exclusive two-pion electroproduction has been studied by ZEUS at HERA in the range 
0.4 < M^^ < 2.5 GeV, 2 < < 80 GeV^ 32 < W < 180 GeV and \t\ < 0.6 GeV^ The 
mass distribution is well described by the pion electromagnetic form factor, [^^^(M^^)^, 
which includes three resonances, p, p'(1450) and p"(1700). 

A dependence of |-F^(M7r7r)P is observed, visible in particular in the interference region 
between p' and p". The electromagnetic pion form factor obtained from the present 
analysis is lower (higher) than that obtained from e~^e~ tt+tt" for Af^^r < 0.8 GeV 
(0.8 < A'U^ < 1.8 GeV). They are equal within errors for M^^ > 1.8 GeV. 

The dependence of the cross-section ratios Rp' = cr(p' — )■ 7r7r)/cr(p) and Rp" = a(p" — )• 
7r7r)/cr(p), has been studied. The ratio Rpf rises strongly with Q^, as expected in QCD- 
inspired models in which the wave-function of the vector meson is calculated within the 
constituent quark model, which allows for nodes in the wave-function to be present. 
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Parameter 


ZEUS 


PDG 


Up (MeV) 


771 ± 2tl 


775.49±0.34 


Tp (MeV) 


155 ± 5 ± 2 


149.1±0.8 




-0.27 ± 0.02 ± 0.02 




Up, (MeV) 


1350 ± 20^^!^ 


1465±25 


Tp' (MeV) 


460 ± 30+^° 


400±60 


7 


0.10 ±0.02^°:°? 




Mp// (MeV) 


1780 ± 20^^^ 


1720±20 


Tp// (MeV) 


310 ± 301^^ 


250±100 


B 


0.41 ±0.03 ±0.07 




n 


1.30±0.06l|J:i^ 





Table 1: Fit parameters obtained using Ft^^Mt^t^) parameterization. Masses and 
widths are in Me V. The first uncertainty is statistical, the second systematic. Also 
shown are the masses and widths from the PDG l4^ . 



g2(GeV2) 


2-5 


5-10 


10-80 




-0.249 ± 0.008+°;°°^ 


-0.282 ± 0.008+°:°°^ 


-0.35 ±0.02 ±0.01 


7 


0.100 ±0.009 ±0.003 


0.098 ± 0.012t°;[!°| 


0.118 ±o.o22j:[!:°[!^ 



Table 2: The dependence of the /3 and 7 parameters. Masses and widths are 
fixed to the values given in Table [21 The first uncertainty is statistical, the second 
systematic. 



Q2(GeV2) 


2-5 


5-10 


10-80 


Rpi 


0.063±0.006± 0.004 


0.081±0.007l[]:[]g^ 


0.122±0.008lJ]:[Jg^ 


Rp" 


0.027±0.006+|]:|J[;^ 


0.026±0.006± 0.003 


0.039±0.010lJ^;[![J^ 



Table 3: The dependence of the ratio Ry for V = p' and p" . The first 
uncertainty is statistical, the second systematic. 
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Figure 1: Comparison between the data and the Zeusvm MC distributions for 
the transverse momentum, p^, of tt"*" and 7r~ particles for different ranges of 
and Mjr,r as indicated in the figure. The MC distributions are normalized to the 
data. 
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Figure 2: Comparison between the data and the Zeusvm MC distributions for 
Q^, W , \t\, cosOh, $/i and (ph for events within mass range 1.1 < M^^r < 1-6 GeV. 
The MC distributions are normalized to the data. 
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Figure 3: Comparison between the data and the Zeusvm MC distributions for 
Q^, W , \t\, cos6h, $/i and (ph for events within mass range 1.6 < M^-tt < 2.1 GeV. 
The MC distributions are normalized to the data. 
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ZEUS 




Figure 4: The two-pion invariant-mass distribution, Mtttt, where A^^^-tt is the 
acceptance- corrected number of events in each bin of 60 MeV. The dots are the data 
and the full line is the result of a fit using the Kuhn-Santamaria parameterization. 
The dashed line is the result of the pion form factor normalized to the data and the 
dash-dotted line denotes the background contribution. 
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Figure 5: The two-pion invariant-mass distribution, M^^j^, where A^^^r is the 
acceptance corrected number of events in each bin of 60 MeV. The dots are the data 
and the full line is the result of a fit using the Kuhn-Santamaria parameterization. 
The contributions of the three resonances p, p' and p" are shown as dashed, dash- 
dotted and dotted lines, respectively. The sum of their interferences is shown by the 
long- dash- dotted line. The background is presented as the sparse dotted line. 
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Figure 6: The acceptance- corrected cos 6^ distribution for different M^^^^ intervals, 
with the mean mass values indicated in the figure. The lines represent fits to the 
data as discussed in the text. 
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Figure 7: The fitted parameter r as a function of the two-pion invariant mass, 
Mtttt- Only statistical uncertainties are shown. 
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Figure 8: The two-pion invariant-mass distribution, Mj^^, where N^^^^ is the 
acceptance- corrected number of events in each bin of 60 MeV, for three regions of 
, as denoted in the figure. The dots are the data and the full line is the result 
of a fit using the Kuhn- Santamaria parameterization. The dashed line is the result 
of the pion form factor normalized to the data and the dash-dotted line denotes the 
background contribution. 
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Figure 9: The pion form factor squared, l-FVp, as a function of the 7r~^7r~ invariant 
mass, Mtttt, as obtained from the reaction e~^e~ — t- tt+tt" /7g[[73| , [7^[7^[7g| /. The 
shaded bands represent the square of the pion form factor and its total uncertainty 
obtained in the present analysis for three ranges of Q^: 2-5 GeV^ (crossed lines), 
5-10 GeV^ (horizontal lines) and 10-80 GeV^ (vertical lines). 
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Figure 10: The pion form factor squared, I-FttP, in the p mass region, as a 
function of the vr+vr" invariant mass, Mj^^, as obtained from the reaction e~^e^ — t- 



fT5 [\jmjmj^\J^ . The shaded bands represent the square of the pion f( 



orm 



factor and its total uncertainty obtained in the present analysis for three ranges 
of Q"^: 2-5 GeV^ (crossed lines), 5-10 GeV^ (horizontal lines) and 10-80 GeV^ 
(vertical lines). 
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Figure 11: The ratio Ry as a function of for V = p' (full circles) and p" 
(open squares). The inner error bars indicate the statistical uncertainty, the outer 
error bars represent the statistical and systematic uncertainty added in quadrature. 
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